Spin current as a response to external stress 



Prashant Sharmg0 
Materials Science Division, Argonne National Laboratory, Argonne IL 60439 
(Dated: February 2, 2008) 

It is theoretically predicted that a traveling shear wave will create a spin current in certain direct-gap 
(for example III-V compound) semiconductors with contributions from both the valence bands and 
the conduction band (for n-doped semiconductors). We show that this spin-current is a property 
of the Fermi-Dirac sea, and is controlled by a geometric phase accumulated by the strain-induced 
Rashba parameters in a cycle. 



In recent years there has been a great deal of exper- 
imental and theoretical interest in the manipulation of 
electronic spin in semiconducting materials Gener- 
ating spin current in a controlled way has emerged as 
one of the central problems of spintronics 0, |^ . Early 
interest in the obtaining spin polarized current by adia- 
batic pumping through quantum dots |^] was not with- 
out experimental difficulties More recently, progress 
has been made in generating spin currents of fixed po- 
larization direction in bulk semiconductors 6] using the 
extrinsic spin-Hall effect H ? ^^^^ two-dimensional 
layers of p-GaAs [loj using the intrinsic spin-Hall ef- 



fect 



rs oi p-l^a. 
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In contrast to charge current, spin angular momen- 
tum current can occur without net transport of parti- 
cles. Therefore one expects to find a method for creating 
a spin current without using external electric fields which 
lead to Fermi surface effects. In this paper we show that 
it is possible to use the coupling of electrons to local de- 
formations in the crystal in such a way as to transfer 
external angular momentum to the electron's spin [2^ . 
This offers an alternative to the aforementioned methods 
for generating and controlling spin currents. 

Under shear stress that breaks crystalline symmetry lo- 
cally, the motion of a Bloch electron is affected depending 
on its spin state. In a semiclassical picture, two different 
effects occur together: (i) Because of local lattice defor- 
mations that break inversion symmetry, if a particle mov- 
ing in (for example) the positive z direction enters such a 
region and acquires a deflection in the positive x-direction 
(dashed line in Fig. 1.), then the oppositely moving par- 
ticle (solid line in Fig. 1) is deflected with a different (in 
magnitude) x component of the velocity, (ii) Because of 
strain-induced spin-orbit interaction the quasiparticle's 
energy in this new environment depends on its spin di- 
rection and is odd in its velocity. Its effect is, therefore, 
the same as a spin-dependent vector potential for the 
quasiparticle. If the strain is position dependent along 
the X-axis (for example) so that the vector potential has 
a non-zero curl, then the resulting orbital magnetic field 
is spin-dependent, but insensitive to the direction of ve- 
locity. Let us say the situation is such that the magnetic 
field has only a ?/-axis component. The Lorentz force 
that acts on quasiparticles moving in the z-direction (see 



Fig. 1.) is such as to deflect the up-spins (blue dashed 
line) in the positive x-direction and the down-spins (red 
dashed line) to the negative x-direction. Under a reversal 
of velocities (solid lines in the figure) an opposite deflec- 
tion will occur. Because of the asymmetrical scattering 
(i), however, these two events will not cancel each other. 
Therefore, as this lattice distortion moves through the 
crystal it will create a net spin angular momentum cur- 
rent flowing in the x-direction. 




FIG. 1: Schematic of time reversed semiclassical trajectories 
of spin up and down particles. The presence of local lattice de- 
formation and the induced spin-orbit interaction breaks time- 
reversal symmetry and creates unequal deflection of these tra- 
jectories leading to a net spin current. 

To put the above arguments on a firm footing, let us 
consider a direct gap semiconductor crystal with cubic 
and time-reversal symmetry. The Hamiltonian for con- 
duction band electrons conserves the spin S and is writ- 
ten in the spherical approximation as 
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where rric denotes the effective mass in the conduction 
band and k is its crystal momentum. The Hamiltonian 
for the valence band is more involved if the bands are 
degenerate. For definiteness, we will consider the case, 
where at the zone center (the F point) the valence band is 
triply degenerate (without spin) and well-separated from 
the conduction band. The valence band Hamiltonian 
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is therefore described by the orbital angular momentum 
projection / = ±1,0 with a (negative) reduced effective 
mass for the electron. Because of spin-orbit coupling 
the conserved quantity is the total angular momentum 
J = L -\- and the helicity k • S/k. Therefore the or- 
bital angular momentum bands are further split into the 
(doubly degenerate) heavy hole (HH) and light hole (LH) 
bands with J = 3/2 and helicity quantum numbers ±3/2 
and ±1/2, and the J = 1/2 and helicity ±1/2 SH-band. 
The Hamiltonian near the k = (F) point is given by 
the Luttinger Hamiltonian [l^ : 
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where m is the free electron mass and the 7s are Lut- 
tinger 's parameters, and the spin S is in the spin-3/2 
representation of SU(2). In the spherical approximation 
the energy spectrum of the various bands is obtained as: 



2m. , 



, /i = HH, LH, SH, c. 
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where denotes the effective mass of the electron in 
the /i-th band. Let us note here that the spherical ap- 
proximation comes with an upper cut-off that we denote 
by kc. 

When the crystal is subject to a transverse (shear) 
wave, both inversion symmetry as well as time-reversal 
symmetry (because of the traveling wave-form) are bro- 
ken. To consider transport along one direction we choose 
an acoustic wave traveling along the x-axis with parti- 
cle motion along the y- and z-axes so as to create time- 
dependent strains Sxy = (l)z{x^t) (along the z crystalline 
axis), Eyz = (j)x{x^t) (along the x-crystalline axis), and 
Sxz = (j)y{x,t) (along the y crystalline axis). To de- 
scribe the response of the electrons to this perturba- 
tion we adopt the semiclassical approach of Ref. in 
which we consider a wave-packet centered at the coordi- 
nate X = (x, y, z) at a given time with its spread small 
compared to the wavelength 27r/g of the acoustic pertur- 
bation. The local Hamiltonian for such a wave-packet 
can be written down in the basis state of the Bloch wave 
|6'/ik, X, t). In the "clean" limit — when the wave- vector 
of the sound wave q > l//p, the inverse particle scatter- 
ing length — the electrons move in phase with the travel- 
ing shear wave. The dominant contribution to spin-orbit 
interaction from such a strain is linear in the crystal mo- 
mentum k and affects both the conduction and valence 
bands so that the energy of the wave-packet is 
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where the spin of the quasiparticle is and the vel ocity 
Cjj, depends on the semiconductor band structure |l5j . 



Besides this Rashba correction, iJ^^, the strain also cre- 
ates a deformation potential [l^ [13, llM • 
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) , and we have written the 



where m,,^ 

electron-phonon deformation potential in the deformable 
ion approximation 23]. With these approximations we 
obtain the energy of the electrons in the /i-th band per- 
turbed by the acoustic wave 



^ H;{t), /i = c,HH,LH,SH (If) 

The time-dependent perturbations H^^^{t) ± H^{t) 
change the crystal momentum kx according to the semi- 
classical Bloch equations of motion: 

dkx 
~dt 
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Note that the wave packet's velocity transverse to the 
wave propagation direction also changes: 



fi/koj 



dy 
dt 

dz hk 
dt m 



h \ Jy rUa 



hkz 



hkr 



h 



fxikqi 



I \ Iz m,. 



Jl/kq' 



^y ~ 



(3a) 
(3b) 



Eqn.[21can be readily integrated to give 

k^t) = 1^ ^{el'' ^*"^(*") , E{t')} + e/o '^*"^(*") 

where kxo denotes the crystal momentum in the bare 
(unperturbed) crystal which is the situation for time t < 
0, and { A, B} denotes the symmetrized product. 

The current response to this time-dependent pertur- 
bation can be calculated analytically if we restrict our- 
selves to terms up to 0^ <C 1. This scheme requires 
that the strains satisfy {vc/vs)\(l)a\ 1, where Vc = 
max[C^, fi/cc/m^], kc the upper cut-off, and Vs the ve- 
locity of sound. Upon expanding the exponential in the 
equation above we determine kx{t) as a function of the 
spin S and the equilibrium momenta kxo^ ky^z' 
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The spin S also acquires additional dynamics from the 
perturbation. This is identified as arising from an effec- 
tive Zeeman field B with components (in units of C^) 

Bx — (pykz — (pzky'f By = (pzkx ~ (j^xkz'i Bz = 4^xky — (j)ykx- 

The perturbation gives rise to a (time-dependent) 
Rashba splitting between the S bands 2\S\C^B = 

2\S\C^^ B1 + + B1. and the change in the equa- 
tions of motion of the spin because of the Hamiltonian 
(Eqn.[id|) is obtained using the SU(2) algebra [S^'.S^] = 
iheabcS'': 

f=C,5xB (5) 

The generalized current density in the a-direction of 
charge (6 = 0), or the 6-th component of spin S is given 

by 

fait) = ^ E ^{S^^,-,t\{^,S\t)}\S^^^.,.,t), 

(6) 

where Q is the unit cell volume, 5*^^^ = e, the electron 
charge, and n^(ko) denotes the quasi-particle occupation 
in the equilibrium state with momentum ko, band index 
/i, and with quantum number and |S'/ik, x, t) denotes 
the Bloch wave at coordinate x and the time t. We have 
not included the S- index for the Fermi occupation num- 
ber since in equilibrium these states are degenerate and 
are equally occupied. Because the momenta ky and kz are 
conserved there are no charge currents along the y- and 
z-directions. The only possibility of a non-zero charge 
current is along the x-direction because kx is not con- 
served. However, the change in kx (see Eq. EJ does not 
create any charge current as E and the explicitly spin- 
dependent part of F sum to zero over the 6'-bands, while 
the product of kxo and the spin-independent part of F 
averages to zero over the zone. Furthermore, there is no 
spin current along the y- and z-directions because the 
momenta ky^z are uncorrelated with spin. We therefore 
focus on the spin current 

m= E ^^{Sf^k,^,t\{^,s\t)}\Sf.k,^,t), 

/i,ko,5 ^ 

(7a) 

that can be calculated using Eq. ^ Only those terms 
in the summand contribute that are: (i) functions of 
/c^Q, ky^ /c^, and (ii) quadratic in the spin S. Therefore 
we look at those terms in Q that are functions of /c^q 
and are linear in S. Only the last term in Eq. ^ con- 
tributes, and we focus now on simplifying it. In doing so 
we shall make use of the smallness of the strain param- 
eters \(j)a\ <C 1: Because {E{t^)^ F{t^^)} is already second 



order in 0^ we will neglect the additional dynamics of the 
spin (Eqn. [5|) in the evaluation of as its time evolution 
will bring in higher powers of (pa- We thus obtain 




- S«{t)A<PyMt')kl + S'{t')A<pJ4t')kl}) (7b) 

where A(j)a = 4^a{t) — (l)a{t')- The terms we have neglected 
involve time derivatives of the spin operator. For conduc- 
tion band spin this of course vanishes. For the valence 
band, however, ((SJ is not the only dynamics of the spin, 
as only helicity is a good quantum number. The form 
of the Luttinger Hamiltonian ()lb|) implies that the time 
derivative of spin will be quadratic in the spin. Together 
with S^{t) present in (|7b|) this will bring in a three spin 
expectation value in the equilibrium state. This vanishes 
upon doing the helicity sum, thus justifying the neglect 
of spin dynamics in determining the spin current. 

Let us write (l)y^z{x,t) = fy^z(l)x{x,t) gy^z^x{x,t), 
where (j)x is out of phase with (j)x, and fy^z^ 9y,z are time 
independent. Using this we obtain the dc spin current: 

£=si:^|p||*'**"'(^-^) («) 

ko,/i * ^ ^ 

•viz) _ n^(ko) h^kl . royjz)'^. x 

where the horizontal bar denotes time average over the 
wave's period. We immediately see that the dc spin 
current is determined by the angular momentum im- 
parted to the quasi-particles along the polarization axis 
by the traveling wave. Within the summand, the quan- 
tity = — 2m~^ . Over most of the zone the 
semiconductor band mass is much smaller than the elec- 
tron mass and we can use ^ —2m~^. It follows 

/J/ 

that the net spin current is independent of the sign of the 
effective mass: All bands contribute together. 

The spin current is controlled by the geometric phase 
that a quasiparticle accumulates (in a period) as it moves 
in a trajectory tracking the ^-parameter. The Rashba 
Hamiltonian (|ld|) implies a spin-dependent (non-abelian) 
vector potential for the fi-th band: A = m^{C^/h)(j) x S. 
The corresponding magnetic field is 

= 0,By= m^{C^/K)[d,<PyS'' - d,<P^S% 
B, = m^{C^lh)[d^<t>.S^ -d^<t>,S']. (9) 

The motion of the particle in the x — z plane is controlled 
by (j)y so that the flux of this magnetic field through a 
unit area oriented in the y direction is proportional to 
By(j)y. Similarily, the motion of the particle in the x — 
y plane is controlled by (j)z and the corresponding flux 
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is proportional to Bz(l)z' Using Q and comparing the 
resulting expression with the spin current ((7b|) it is seen 
that the difference of these two fluxes determines the spin 
current that is proportional to ii[S^{By(j)y — Bz(l)z)]' 

The magnitude of the spin current can be estimated 
as follows. Consider the case when (j)^ = and let us 
also choose fy = = gz and fz=gyy^O,so that the dc 
current is only polarized in the z-direction: 



ii{S'^)/h (10) 
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Such a choice of shear strains can be obtained by ap- 
plying shear stresses schematically drawn in Fig. 2. The 
contribution of the integral in ([Tn|) comes from parts of 
the zone where the effective mass is significantly different 
from the electron mass. To estimate a lower bound for 
the current let us limit the summation range to an up- 
per cut-off wavevector /cc, alluded to earlier, and consider 
n-doped semiconductor at low temperature. The lower 
bound is estimated considering only the conduction band 
(/i = c) for n-doped semiconductor: 



that is essential is that, besides crystal momentum, the 
electron's total angular momentum J is conserved. The 
perturbation alters the crystal momentum along one (x) 
direction thereby leaving only J in that direction con- 
served. Thus dJx/dt = = h{v X /c)aj + dSx/dt. Using 
Eqns. ini and ^ the first term, which is proportional to 
the spin current, can be calculated and its projection onto 
S^{t) gives Eqn.[2b|for spin current. 

In conclusion let us remark that unlike the intrinsic 
spin-Hall effect 0, 0, ^ spin-acoustic effect de- 
scribed here arises from explicitly breaking time-reversal 
symmetry. This leads to contributions from all occupied 
bands that are affected by strain-induced spin-orbit in- 
teractions. In this sense the spin current can be thought 
of as a diamagnetic response of the semiconductor to ef- 
fective spin-dependent magnetic fields created by non- 
uniform stress. 
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using GaAs data: Cc = 8 x lO^cm/s mass ratio 

rric/m = 0.05 that is fairly representative of most semi- 
conductors, and taking /cc = 3 x lO^cm"^ = 10~^A~^ 
{kc = kp for low doping with Ep = 3 meV, Ac is the elec- 
tron's Compton wavelength), and the velocity of sound 
Vs = lO^cm/s. The value of strain (0^)^^^ will depend on 
how the acoustic wave is generated. One way to generate 
such a wave is to contact the crystal by a piezoelectric 
material in which acoustic waves are induced by applied 
ac electric field. This conventional technique gives (for 
electric field strength of lOV/mm and using GaAs elastic 
and piezoelectric constants) a strain (^^)^/^ = 10~^ [l^ . 
Furthermore, we are assuming that the shear wave has a 
wavelength comparable to, or smaller than the momen- 
tum scattering length Ip. At low (liquid nitrogen) tem- 
peratures, the value for Ip ^ Ifim in semiconductors at 
low doping . Acoustic waves of such wavelengths have 
been obtained in semiconductor crystals [i^. We also 
note that the spin lifetime Ts ^ ns ^lla in several mate- 
rials, which allows experiments at GHz frequency. Using 
these parameters we obtain {e/h)j^ ^ 10^ Amp/cm^, 
which is not unreasonable, especially since a spin cur- 
rent smaller than this magnitude has been detected 0, 
and also because the contribution from the valence band 
comes from the entire zone and is therefore much larger 
than this lower bound. 

Our calculation above may seem to be limited by the 
choice of an unperturbed Hamiltonian that has spheri- 
cal dispersion in the valence band. It should be empha- 
sized, therefore, that the only feature of the valence band 
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